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The Motions of the Stars in a Kapteyn- Universe. 

By J. H. Jeans. 

1. In his recent paper “ On the Distribution of the Stars in Space,’ 5 * 
and again in his “ First Attempt at a Theory of the Arrangement and 
Motion of the Sidereal System,”! Professor Kapteyn has given a 
masterly summary of our present knowledge, mainly due to his own 
researches, of the arrangement of the stars in space. In brief, he 
finds that the arrangement of the stars may be regarded, with very 
tolerable accuracy, as consisting of a system of spheroidal shells, each 
of constant density. The spheroids are similar and similarly situated; 
they are oblate, having their axes in the ratio 5*102 : 5*102 : 1, their 
equal longer axes being in the galactic plane and their common centre 
being not far removed from the position of the sun. 

In the second of these papers Professor Kapteyn attempts a 
dynamical investigation of such a system of stars. He assumes that 
the stars are in a dynamically steady state, and makes the further 
supposition that along the shortest axes of these spheroids ( i.e. in the 
directions of the galactic poles) the law of distribution is Maxwell’s 
well-known law for the distribution of molecular velocities of a gas in 
isothermal equilibrium. 

Kapteyn’s first assumption can neither be proved nor disproved * 
in the present state of our knowledge it is probably the most reasonable 
assumption to take as the basis of a dynamical discussion. But, for 

* Astrophys . Journ 52(1920), 23. + Ibid. (1922). 
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reasons which will appear later, the second assumption seems much 
less plausible. Indeed, in introducing it Kapteyn writes : ({ We will now 
introduce some considerations from the kinetic theory of gases, the 
applicability of which to the stellar system might be considered 
doubtful. . . . The results which will be derived cannot lay claim to 
be demonstrably correct.” 

The main object of the present paper is to show that the introduction 
of Kapteyn’s dubious second assumption is unnecessary. As we shall 
now see, the principal features of the motion of the stars can be deduced 
from a knowledge of their distribution in space, combined only with 
the single assumption that the system of stars is in a steady state. 

2. Kapteyn’s given distribution is one of symmetry about a plane 
(the Milky Way) and about an axis (the pole of the Milky Way). 
Taking the centre of the system as origin and the plane of the Milky 
Way as the plane 2 = o, we introduce cylindrical co-ordinates ST, $, z , and 
denote the velocities of a star in the directions of ST, 0 , z increasing by 

n,©,z. 

Let the stars be of various distinct types, distinguished by the 
suffixes 1, 2, 3. ... 

Let the law of distribution of stars of type 1 at any instant t be 
assumed to be 

/^n, ©, Z, ST, 0, ?, £)dll d © dZ dTS ST dOdz . . ( 1 ) 

this representing the number of. stars in the element of volume 
dfcr ST dO dz of which the velocity components at time t lie within the 
small range between II and II + dll, © and © + d©, Z and Z + dZ. 

Each star will be in motion under gravitational forces which may 
be regarded as made up of two parts— 

(A) The forces from adjacent stars. 

(B) The forces from the gravitational field of the universe as a 
whole. 

A simple calculation shows that in our own system of stars forces 
(A) are in general quite insignificant in comparison with forces (B).* 
Thus we may expect a system of stars such as our own to tend first to 
assume a state of steady motion under forces (B). When such a state 
has been reached the average effect of forces (B) becomes zero, and 
forces (A) now come into play, until gradually the system assumes a 
final steady state under both forces (A) and (B). This final steady 
state is exactly analogous to Maxwell’s steady state in the kinetic 
theory of gases; under it the velocities of the stars in any small region 
of space conform to Maxwell’s law of distribution of velocities, so 
that these can show no preference for special directions in space, and 
there can be no star-streaming. The fact that star-streaming is observed 
in our universe affords conclusive proof that this final state has not yet 
been attained, a conclusion which is fully in accordance with calculations 
as to the time which would be required for its attainment.! 

On the other hand, the time required for the attainment of a steady 
state under forces (B) alone is of an altogether lower order of magnitude, 
so that it is quite possible that a steady state of this kind has already 

* M.N., 74(1913), 109. f Lx. ante , p. 112. 
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been attained. In agreement with Kapteyn, we shall assume this to be 
the case throughout the remainder of the present paper. 

3. For such a steady state it is readily shown * that the law of 
distribution (1) assumes the special form 

/ X [J(n 2 + © 2 + Z 2 ) - Y, tt®]dll d © dZ dV5 73 dOdz . ( 2 ) 

where Y denotes the gravitational potential of the universe as a whole. 

In this formula, f x is no longer specified as depending explicitly upon 
t , but it will be understood that the action of the feeble forces (A) 
results in a slow secular change in the constants which occur in the 
specification of / r It is easily shown that this formula demands star¬ 
streaming in the direction of ©, the direction which Charlier and 
Kapteyn assign as the axis of star-streaming in our system. 

Let Vj denote the number per unit volume of stars of type 1 in any 
small region of space, and let (II©)! denote the average value of the 
product II© for these stars. Then, under the law of distribution (2), 

^(11®), = j j JYn© cm d® dZ = o, 

and similarly 

■' 1 (©Z ) 1 = o, > , 1 (zn ) 1 = o, 

Let us introduce quantities p v q 1 defined by 

Vl (a^v(Z\= Pl , . . . (3) 

Taking u, v, w for the moment to represent components of stellar 
velocity along any three rectangular axes x, y , z, it is easy to deduce- 
from the law of distribution (2) three equations of the type f 

jhWi]+^( 4 ] 4 WH]= • • (4) 


these being, in point of fact, the hydrodynamical equations of motion for 
the v x stars of the first type in a unit volume at x, y , 3, taking into 
account the fact that forces (A) are negligible. 

Transforming to the cylindrical co-ordinates trr, 0 , z already employed, 
these equations become 


dz 1 dz 

• 

• 

• ( 5 ) 

dq, dV 

dO~ Vl dO * 

- 

• 

• ( 6 ) 

9 /r _. \ ^ 0 Y 

• 

. 

• ( 7 ) 


Since the field of stars is assumed to be symmetrical about the axis 
of z, 0 Y /30 must vanish, and equation (6) merely expresses the obvious 

* Problems of Cosmogony and Stellar Dynamics } p. 236. 
t Ibid. } p. 232. 
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fact that q 1 cannot vary with 0 . It is convenient to rewrite equation (7) 
in the form 




Pi-to 



( 8 ) 


and equations (5) and (8) now contain the material for discussion. 

The equations, as at present written, refer only to stars of one special 
type, hut we are free to add together corresponding sides of any number 
of such equations. In this way we obtain equations such as 


dp dV 

- = v — 
dz dz 


(9) 


d P . p-q^ dY 


(10) 


where p, q , v now refer to all stars of any type or class desired, e.g. all 
B-type stars, or all luminous stars or all dark stars. But Y is in every 
case the gravitational potential of all stars and of all dark matter in the 
universe. It will be noticed that equations (9) and (10) are no longer 
hydrodynamical equations. 

Kapteyn’s investigation gives us in effect the values of v at every 
point of the universe. If we assume the stars to have everywhere the 
same known average mass, it is a mere matter of calculation to determine 
Y at every point. This being done, equation (9) enables us to find p at 
each point, and equation (10) then determines q at each point. Since 
v is already known, a knowledge of p and q gives the mean values of 
II 2 , Z 2 , and © 2 at each point of space and for any special class of stars we 
please. The mean values of II 2 and Z 2 will give us information as to 
the mean peculiar motions of the stars, while the mean value of © 2 - Z 2 
will give information as to the extent of star-streaming. 

4. As a preliminary illustration of the use of our formulae let us 
consider their application to the B stars. Charlier * has shown that 
these are confined to a flat, biscuit-shaped space surrounding the galactic 
plane. The sun is only about 90 parsecs from the centre of this 
region, its radius in the plane of the galaxy being about 500 parsecs 
and its extension on either side of this plane only about 150 parsecs. 
As a matter of convenience let us suppose the boundary of this region 
to coincide with the first ellipsoidal shell of Kapteyn, of which the 
boundary is 


6X 2 ?? _ 

a ~ 2 + ^~ 


h 


1 ( a = 602 parsecs 

where < 0 

\ c = 118 parsecs 


00 


Let us further, as an approximation which appears to be in very fair 
agreement with the facts, assume the B-type stars to be arranged in 
similar ellipsoidal shells inside this one, so that the density v of B-type 
stars is of the form 




• ( i2 ) 


* Lund Meddel 2 (1916), No. 14. 
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Inside the ellipsoid occupied by the B-type stars the potential of the 
whole universe, as arranged by Kapteyn, is of the form 

V = C - Aar 2 — Bz 2 , 

where the last two terms arise from the potential of the ellipsoid itself 
and all external ellipsoids add only a constant term to Y. 

From equation (9) we now have 

p -l v % d *=-K i - S-5) 2Bztfe • • (i3) 

If we put Jc)>(x) dx = \p(x) and adjust the additive constant so that 
1 j/(x) vanishes when x = o, then equation (13) becomes 




■ ( 14 ) 


no function of G 7 being added, since p must vanish at the boundary (n). 
Equation (10) now gives 


= ?F(A a 2 - B c 2 )J 1 - ^ 

ZS a 2 v / \ a 2 c 2 / 


or 


< Lz£. = 2(Aa 2 - B6 2 )^ 


■ ( 15 ) 


The maximum value of p occurs at the centre of the system, being 
Bc 2 i/r(i) or Bc 2 j dx. If we assume, for illustration, a law of density 






then <f>(x) = v 0 (x) and 1/^(1) = Jv 0 . In this case the value of pjv at the 
centre is ^Bc 2 . Had we assumed uniform density we should have had 
pjv = Be 2 at the centre. 

In his numerical calculations Kapteyn takes the parsec as unit of 
length and the gravitational mass of the “ average star ” as unit of mass. 
In these units he calculates the gravitational forces at the points (a, o) 
and (o, c), and finds 

2k.a— 33*19, 2Bc = 40*o6. 

In these units the value of p/v at the origin is 

p/v = ^Bc 2 = 1182 . 

To reduce to ordinary C.G.S. units we have to multiply by yM/Z. Here 
l is the number of cms. in a parsec (3'o8 x io 18 ), y is the gravitation 
constant (6-66 x io~ 8 ), and M is the mass of the u average star” in 
grammes. More exactly this last quantity should be the mass of the 
average visible star increased by the ratio of dark matter to visible stars 
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in the universe. Later on we shall determine the value of M as 
6*34 x io 33 grammes. We shall use this value here, so that 

-*y- = 1*37 x io 8 , 

and 

^ = £Bc 2 = 16*2 x io 10 (cms. per sec.) 2 . 

This quantity is of course equal to the mean value of Z 2 or II 2 for all 
the B-type stars near the origin. We see that even close to the origin, 
where they attain their maximum values, Z and II are only of the order 
of 4 kms. a second (or 6 kms. a second in the extreme case of uniform 
distribution). It is at once clear that the low velocities of the B-type 
stars are fully accounted for by their observed distribution in space: it 
is not necessary to interpret them as a consequence of their supposed 
exceptional mass. 

Independently of the law of distribution of density, equation (15) 
gives 

q— p ftf \ 2 10/GT\ 2 . 

±—±-= 15,252^—^) =209 x.io iU l —j (cms. per sec.)* 5 . 

This quantity represents the mean value of © 2 — Z 2 or © 2 — II 2 . When 
it has a positive value it is clear that there must be star-streaming in 
the directions of 0 increasing and decreasing. If IT, U" are the velocities 

of the two star streams, we shall see later (§ 6) that 2 —E is equal to 

IT U". Thus in the special case now under discussion the velocities of 
the star streams are of the order of 14*5 kms. a second at the boundary 
£7 — a (602 parsecs) and decrease uniformly with ZS, precisely as though 
the two star streams were rigid bodies rotating in opposite directions. 

At Charlier’s position for the sun (GT = 90 parsecs), the velocity of star¬ 
streaming for the B-type stars will be only about 2*2 kms. a second. 

5. Let us now apply the same method to the general mass of the 
stars. From equation (9) we have 

p =~Fi dz - • ■ • • (l6) 

where the integration is along a line parallel to the axis of 2 until a 
region is reached at which the star density v is negligible. The integral 
can only be evaluated by quadratures. Kapteyn’s investigation gives 
the value of v at every point, and he has calculated the value of — dY/dz, 
which he denotes by G cos </>, at certain selected points. I have first 
evaluated p at points along the galactic axis (p = o), taking the thick¬ 
nesses of Kapteyn’s various “ shells ” as the elements of integration in 
equation (16). The calculation is exhibited in Table I. The entries 
Ex. 1, Ex. 2, etc., refer to shells obtained by extrapolation beyond the 
ten shells of Kapteyn. No estimate is made of the density or thickness 
of these shells, the extrapolation being only with reference to their con¬ 
tribution to p. The last column gives the square root of the mean value 
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of Z 2 or II 2 in kilometres per second on the assumption that M has the value 
already assigned to it in § 4. This may he regarded as giving approxi¬ 
mately one component, other thai; ©, of the star’s peculiar velocity; it is 
seen to range from about 13 kms. a second in the close neighbourhood 
of the origin, through a value of 17 kms. a second at about 300 parsecs 
from the origin, to a final value of about 23 kms. a second, which is 
approximately attained at about 1600 parsecs from the origin. 


Table I. 

Values of p and Z along Galactic Pole. 


Shell. 

Mean v. 

0V 

Mean • 

dz in 
parsecs. 

0V 

v^dz. Ellipsoid. 

z in 

parsecs. 

p- 

V. 

p 

z 

(kms. per se 

I 

•036 

22 ’0 

118 

93‘5 

Centre 

0 

504-0 

•0451 

11,200 

12*4 

2 

•023 

48'0 

80 

88-3 

I. 

118 

4 i °‘5 

‘0285 

14,400 

14-0 

3 

•014 

58*5 

98 

803 

II. 

198 

322*2 

*0180 

17,900 

157 

4 

•0090 

64*8 

n 7 

68-2 

III. 

296 

241*9 

•0113 

21,400 

17*1 

5 

•0057 

66*3 

140 

52*9 

IV. 

413 

1737 

•00715 

24,300 

18*2 

6 

•0036 

65*0 

164 

38 ’4 

V. 

553 

120*8 

•00451 

26,800 

19-2 

7 

•0023 

62 3 

185 

267 

VI. 

7 i 7 

82-4 

•00285 

29,000 

19*9 

8 

•0014 

5 »-s 

212 

17*4 

VII. 

902 

55-9 

•00180 

31,100 

207 

9 

•0009 

S4'o 

251 

12 ’2 

VIII. 

1114 

38*5 

*00113 

34,100 

21 *6 

10 

•0006 

487 

295 

8*6 

IX. 

1365 

26-3 

•000715 

36,800 

22'4 

Ex. i 




61 

X. 

1660 * 

177 

•000451 

39,3oo 

23 2 

Ex. 2 




4*2 

Ex. I. 


II-6 

•000285 

40,700 

23-6 

Ex. 3 




2 *9 

Ex. II. 


7*4 

•000180 

41,100 

237 






Ex. III. 


4*5 

•000113 

40,000 

23*4 


The next stage in 

the 

computations was 

to construct tables giving 


the values of p throughout the whole system. The method adopted 
was to regard p at any point, as before, as arising from contributions 
from all the shells standing outside this point. It was found that the 

0 y 

contribution v—dz from any specified shell did not vary greatly for 

different points of this shell until very low galactic latitudes were 
reached. For instance, for shell 2, if the contribution in latitude 90 
is taken to be i*ooo, those in latitudes 75 0 , 6o°, 45°, 30°, 15°, and o° 
are respectively *999, *996, *990, *97, *90, and -67; for shell 8 the 
corresponding numbers are ‘984, *950, -86, 745, -57, and *27. From 
a table giving the ratio of contribution in any latitude to that in 
latitude 90° it was found easy to evaluate p at any desired point. The 
following table gives the values found fory? and pfv at a series of points 
in the galactic plane, the points selected being the boundaries of the 
various shells. The table also contains the value of Z, the square root 
of (Z 2 ) in kilometres per second, the value of (q - p)jv calculated 
from equation (1 o), and the square root of (® 2 - Z 2 ) in kilometres 
per second. 
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Table II. 

Values of p, Z, q , and 0 in Galactic Plane. 


Ellipsoid. 

XX in 
parsecs. 

V. 

P- 

P' 

V 

q-p 

V 

Z 

(kms. per sec.). 

V(© 2 -z 2 ) 

(kms. per sec.). 

Centre 

O 

•0451 

504 

11,200 

0 

124 

O 

I. 

602 

•0285 

368 

12,900 

14,300 

13 *3 

14*0 

II. 

1010 

*0l80 

251 

13,900 

25,900 

13*8 

18*8 

HI. 

1510 

•OII3 

165 

14,600 

38,500 

14*1 

23*0 

IV. 

2106 

•00715 

102 

14,400 

51,500 

14 ’O 

26-6 

V. 

2820 

•OO45 

63 

14,000 

62,300 

13-8 

29*2 

VI. 

3656 

•OO285 

36 

12,600 

69,100 

13 1 

30-8 

VII. 

4600 

•.00180 

21 

II,600 

78,200 

12*6 

327 

VIII. 

5675 

•OOI 1 3 

12 

10,600 

83,400 

12*1 

33‘8 

IX. 

6960 

*000715 

7 

9,900 

86,300 

117 

34*4 

X. 

8465 

•OOO45 

4 

8,900 

86,300 

11*0 

34*4 

The tables have given the 

values of p, Z, 

g, © only in 

the galactic 


plane and along the galactic axis. It is easy to obtain a rough idea of 
the values of these quantities in other parts of space. As regards p, 
we have already noticed that the contribution of any individual shell 
to p remains approximately constant until low galactic latitudes are 
reached, so that we may think of p as remaining constant over any 
shell except in low galactic latitudes. In other words, except in these 
low latitudes p and Z may be regarded as functions of v only, and the 
last four columns of Table I. may be regarded as tabulating p and Z in 
terms of v. In low galactic latitudes Table II. may, of course, be used 
in a similar manner. As regards g, we obtain from equations (9) and 
(10) on elimination of Y 

d/q-p\_dp 3 /T\ dp 3 /A 
3 z\ vTZ ) dTZdz\v) dz 3 Xd\v) 

We have seen that in fairly high galactic latitudes p may be 
regarded as a function of v, in which case the right hand of this 
equation vanishes. Again, in low latitudes dp/dz and 3 vjdz are very 
small, so that again the right-hand member is small. We may get a 
rough idea of the extent of star-streaming by neglecting the right-hand 
member of this equation altogether, in which case we find that (g -p)jv 
becomes a function of XS only, and Table II. may be regarded as giving 
this quantity in terms of XS. 

6. If we regard the stars in any small region of space as forming 
two star drifts in the direction of 6 increasing and decreasing, if v, v" 
are the densities of these streams, U', - TJ" the stream velocities, and if 
u', v\ w' and u", v", tv" are the peculiar velocities of individual stars in 
these streams, then 

p — ^^v ' 2 + in unit volume 

— w 2 

q =2 ( u '+ w ') 2 +2 ■(- u "+ u 'y =1/172++ 
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In these equations we have v + v" = v, v'U' = v"U" and u 2 = t/ 2 , so that 
q-p = v'U ' 2 + v"U " 2 ^ wv „ 

V V 

If k v h 2 are the semiaxes of the Schwarzschild velocity ellipsoid at 
any point, we have 

h(f _ 2® 2 q 

p' 

7. In the neighbourhood of our suu we may take approximately 
K/ h 2 = °' 577 , giving V/V = 3 '°°. or 

q — p p 

±—=- = £ X 2 * 00 . 
v v 

In Table II. this relation is seen to be satisfied at G? = 1090 parsecs, 
at which p/v = 14,000 and (q -p)/v = 28,000. On the two-drift hypo¬ 
thesis, the numbers of stars in the two drifts are found to be almost 
exactly in the ratio 3:2, and the relative velocity of the two drifts is 
40 kms. a second. Thus U'=i6 kms. a second and U" = 24 kms. a 
second, giving 

U'U" = 384 x 10™ C.G-.S. units. 

We may notice in passing that the value of U'U" is not very sensitive 
to the ratio in which the stars are supposed divided between the two 
streams. If we had supposed the streams to be equal, we should have 
had U' = U" —20 kms. a second, and the value of U'U" would have 
altered only to 400 x 1 o 10 . 

We have now found that the value of (q—p)jv in Kapteyn units is 
28,000 and in C.G.S. units is 384 x io 10 . Hence, as in § 4, we must 
have 

384 x io 10 = ^ x 28,000, 

V 

whence, on putting y=6'66 x io -8 , 1 = 3^08 x io 18 , we find 

M= 6*34 x io 33 grammes, 

which is the value which has been used throughout the paper, and has 
already been inserted in the values of © and Z given in Tables I. and II. 

This value, being 3*2 times the mass of the sun, is obviously too 
large to represent the mass of a single star. The average mass of binary 
systems is found to be about r6 times the sun’s mass, but these are 
mainly binaries of long period whose origin must probably be traced 
to separate nuclei in the original nebulae, and whose combined mass 
may therefore be expected to be twice that of the average single star. 
Assuming the average star to have a mass of o*8 times the sun’s mass, 
our value of M must be interpreted as showing that if the stars are in a 
steady state, there must be about three dark stars in the universe to 
every bright star. 

8. Subject to our being able to assume the existence of this very 
appreciable amount of dark matter, it appears that the observed motions 
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of the stars can be explained if, and only if, we are free to assume our 
sun to be at a distance of 1090 parsecs from the centre of the galaxy. 
This figure of 1090 parsecs cannot be regarded as having any great 
degree of accuracy, for we have made use of Kapteyn’s observational 
distribution of density, which was derived only on the supposition that 
the sun is approximately at the centre of the system. The true value 
of the sun’s distance from the centre can only be obtained as the result 
of a more refined investigation of stellar density, of which the need is 
fully recognised by Professor Kapteyn, in which full account will be 
taken from the outset of the excentric position of the sun. 

9. Without carrying out this investigation in detail, it is easy to 
discover the general nature of the change to be expected in our estimate 
of 1090 parsecs. In place of the quantities 2, tar, v,£>, q f M so far used, 
let us transform to new quantities 2', £ 7 ', v\ p, q, M' connected with the 
foregoing by the relations 

*2' = k 1 z 9 £7' = ^£ 7 , v — kp\ = q' = k z q } M' = & 4 M 

where k v & 2 , k s , k 4 are constants to which a value will shortly be assigned. 
Transformed to these quantities, our fundamental equations (9) and (10) 
become 

dp _ hje 3 , 0 V' dp p -q _ kjt 3 , 0 V' 
dz~Jc 2 % v dz’’ 0CT' + ©' ~k 2 i k i v to' - 

Thus the fundamental equations transform back to themselves if 
kjc z = & 2 2 & 4 . If this relation is satisfied between the four &’s there will 
be a solution in 2', £ 7 ', etc., corresponding exactly to the solution we 
have already found in 2, £ 7 , etc. To what extent is it possible that this 
new solution should represent the observed facts after allowing for the 
excentric position of the sun? 

The new solution gives a distribution of mass which is of course 
precisely similar to the old, except for proportional changes in length, 
mass, and density. To account for the observed ratio of axes of the 
Schwarzschild ellipsoid, our sun must be at the point in the new 
solution corresponding to £7=1090 parsecs in the old solution; this 
will be given by £ 7 '=1090. To account for the observed density of 
stars near our sun we must have v = 0*0451 in the new solution as 
against v = o*i68 in the old, so that & 2 = 2*68. To account for the 
observed actual velocity of the stars in kilometres per second we must 
have^'M'/v = pM/v, so that we must have & 3 & 4 = ft 2 . The new solution 
would give outer shells of the same shape as those of the old solution, 
but to make the observed density accord with observation a suitable 
value of k x must be chosen. Thus at £ 7 ' = 6960 the new solution gives 
v = *00071 5 x 2*68 = *001910, which observation shows to correspond to 
a distance of £7 = 4500. Thus we must take ^ = 6960-^4500=1*55, 
and it will be found that this value of k x together with the previously 
found value of & 2 = 2*68 gives a new solution which represents the 
observed distribution of stars in the outer layers with tolerable accuracy. 
Eliminating k z from the two relations kjt z = lc^k A and k z k 4 = already 
found, we obtain k 2 = kjk 2 , so that k 4 = 0*76. 

Thus in the new solution the sun’s position (£ 7 ' — 1090 parsecs) 
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will be £7=700 parsecs, and the average mass M will be 4*8 x io 33 
grammes. Taking the average mass of a single star as o*8 times the 
sun’s mass, this indicates two dark stars to each bright star. In the 
new solution the stellar density at the sun’s position (£7 = 700 parsecs) 
will be 0*0451, the value assumed by Kapteyn, and the motion, as 
regards both star-streaming and peculiar motions of the stars, will be 
that actually observed in the neighbourhood of the sun. The new 
solution will represent Kapteyn’s observational material well enough 
as regards the outer shells of stars, but the agreement with observation 
as regards the inner shells needs careful discussion. For reasons which 
will be found in his own paper, Kapteyn concludes that the greatest 
value of £7 which can be reconciled with observation is in the neigh¬ 
bourhood of 700 parsecs. Thus our estimate of £7 is just about at the 
limit permitted by observation. If the estimate had fallen entirely 
outside the limit, we should of course have had to conclude that the 
•stars had not yet attained a steady state of the type assumed in this 
paper. 

In any case we may notice that our estimate of £7 = 700 parsecs 
places the centre of the galaxy right outside the system of the B stars 
as determined by Charlier. This leads us to regard the B stars rather 
as a moving star-cluster than as an integral part of the galactic system 
in steady motion with the rest of this system. We shall return to a 
consideration of this point of view in the next paper. 


The Dynamics of Moving Clusters. 

By J. H. Jeans. 

1. In a recent “ Research on Moving Clusters,” * N. H. Rasmuson 
has provided material for testing a suggestion I put forward some 
years ago t that a moving cluster ought to assume a flattened shape 
as a consequence of encounters with other showers of stars. In 
particular, if the stars met by the cluster are moving at random the 
cluster ought to experience a flattening in the direction of its own 
motion. Certain of the clusters studied by Rasmuson (the Ursa Major 
cluster and that in Perseus, to quote his own instances) show a 
flattening of the type predicted, but others ( e.g . the Scorpio-Centaurus 
cluster) show exactly the opposite effect. 

In my investigation the forces acting on a single star of a moving 
cluster were supposed divided into three parts : 

(a) The forces acting from the remaining stars of the cluster; 

(b) The forces from stars which are undergoing an “ encounter ” 

with this cluster, in the sense of the forces between single 
stars being of appreciable amount; 

( c) The forces from all the remaining stars of the universe. 

* Lund Meddel ., 2, No. 26. t M.N., 76 (1916), 552. 
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